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We generalize the recently developped "internal" Density Functional Theory (DFT) and Kohn- 
Sham scheme to multicomponent systems. We obtain a general formalism, applicable for the de- 
scription of multicomponent self-bound systems (such as molecular systems where the nuclei are 
treated explicitly, atomic nuclei and mixtures of 3 He and 4 He droplets), where the fundamental 
translational symmetry has been treated correctly. The main difference with traditional DFT is 
the explicit inclusion of center-of-mass correlations in the functional. A large part of the paper is 
dedicated to the application to molecular systems, which permits us to clarify the approximations 
that underly traditional DFT. 
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I. INTRODUCTION 

Density Functional Theory (DFT) [lP3 is a tool widely 
used in condensed-matter physics and quantum chem- 
istry to calculate properties of many-electrons systems. 
It is based on the simple local density instead of the less 
tractable iV-body wave function. One of thepillars of 
DFT is the Hohenberg-Kohn (HK) theorem jjj, which, 
in its original form, proves that for any non-degenerate 
system of Fermions or Bosons |l| put into a local ex- 
ternal potential v ex t(r), there exists a unique functional 
of the local one-body density p(r) that gives the exact 
ground-state energy when p(r) corresponds to the exact 
ground-state density. The Kohn-Sham (KS) scheme pro- 
vides a practical method to compute self-consistently the 
ground-state density in a quantum framework, defining 
the local single-particle potential (i.e. an auxiliary non- 
interacting system) which reproduces the exact ground- 
state density @. 

Traditional DFT is particularly well suited to describe 
the electrons in a molecule, but not for the description 
of self-bound systems, such as atomic nuclei, He droplets 
or molecular systems where the nuclei are treated ex- 
plicitly. Indeed, in those systems, external fields are not 
required to obtain bound states. As a consequence, iso- 
lated self-bound systems are plagued by a center-of-mass 
(cm.) problem: for any stationary state, the cm. is de- 
localized in the whole space (because of the translational 
invariance of the Hamiltonian) and the laboratory wave 
function is consequently not normalizable. This prevents 
us from using traditional DFT methods, formulated in 
terms of the laboratory density, when v ext (r) = ||. 
Indeed, the laboratory density is then an indeterminate 
constant |6H8(, which prevents us to construct from it 
a universal functional. Moreover, it is internal proper- 
ties (measured in the cm. frame) and not laboratory 
properties that are of experimental interest. Experimen- 
talists always deduce internal properties using the cm. 



observables (position, momentum or kinetic energy of the 
cm.) and Galilean invariance to transform all the other 
observables into the cm. frame. 

It is thus a question of interest to formulate a rigorous 
DFT formalism and KS scheme in terms of the internal 
density pi„t, having achieved the correct separation of 
internal properties from the cm. motion. Efforts in that 
direction have been made recently in Refs. 0, |9(, but the 
question of a rigorous "internal" KS scheme remained 
open (see the corresponding discussion of Ref. [6| for more 
details). The different approach found in Ref. [Hj, L!J| 
results in a rigorous KS scheme, but that is not "internal" 
(i.e. formulated in the cm. frame of reference), so that 
it is not directly comparable to self-consistent mean-field 
like calculations with effective interactions formulated in 
the cm. frame of reference. 

In Ref. jy], it was proposed to use Jacobi coordinates 
to decouple the cm. properties from internal ones in self- 
bound systems, which permits one to separate the non- 
normalizable part of the wave function from the normal- 
izable part which describes internal properties. More- 
over, an arbitrary translationally invariant potential of 

the form J2i=i v int(ri - R), where R = jjY,f=i r j is 
the total cm. of the particles, was added to the Hamil- 
tonian of a self-bound system. This potential is an "in- 
ternal" potential, i.e. it acts in the cm. frame, and is 
the only form that satisfies all the key formal properties 
@. Of course Vi nt should be zero in the stationary iso- 
lated self-bound case. This is why in Ref. @ the authors 
presented it as a mathematical "auxiliary" to reach the 
desired goal and showed that it can be dropped properly 
at the end, preserving all the conclusions (because the 
internal ground state should by definition remain bound 
at this limit). Nevertheless, its form is suitable to model 
internal effects of fields used in experiments. Through it 
and using Jacobi coordinates it was shown, by a different 
way than those found in Refs. @, Q, the stationary "in- 
ternal" DFT theorem for identical particles: the internal 
many-body ground state can be written as a functional 
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of the internal density pi n t- Then, the corresponding 
"internal" KS scheme was formulated rigorously in the 
cm. frame. This work provided a first step towards a 
fundamental justification for the use of internal density 
functionals for stationary mean-field like calculations of 
nuclei [12| or He droplets [13| with effective interactions, 
showing that there exists an ultimate functional which 
permits one to reproduce the exact internal density. The 
major difference with traditional DFT is that the cm. 
correlations (which appear due to the redundant coordi- 
nate problem) are explicitly included in the functional. 

The aim of the present article is to generalize the inter- 
nal DFT formalism and KS scheme to self-bound systems 
composed of different kinds of particles (i.e. multicom- 
ponent systems) . This is crucial to obtain a fundamental 
justification to the use of density functionals for the de- 
scription of: 

(i) Molecular systems where the nuclei are treated ex- 
plicitly and quantum mechanically. 
Taking into account the quantum nature of the nu- 
clei is important to describe small molecules [14[ 
and solid hydrogen [f5( for instance, and explicit 
treatment of the nuclei is necessary to decribe non- 
adiabatic phenomena. Moreover, taking into ac- 
count the electrons-nuclei correlations can be im- 
portant for the description of some physical proper- 
ties |16j |. For those reasons, various multicomponent 
DFT formalisms have been developped @, Il6l - l20| . 
But in none of them has the cm. motion been sepa- 
rated properly from the beginning, leading to formal 
difficulties (that will be detailed in £jIIID[) . Here, we 
propose a general formulation that overcomes those 
difficulties and is suited for the description of all 
self-bound systems. We show that the application 
to molecular systems is particularly interesting for 
understanding the approximations inherent to tra- 
ditional DFT and thus how to eventually improve 
its results for molecular systems. 

(ii) Atomic nuclei. 

The application of multicomponent internal DFT to 
protons and neutrons gives a further step towards 
the use of mean-field like calculations with effective 
interactions, as currently done in nuclear physics 
|12l | . and shows that the cm. correlations can be 
included in the functional. It thus opens the way to 
an alternative to the numerically very costly pro- 
jection techniques used in nuclear physics to restore 
Galilean invariance [12, l2l| - |23j : they wouldn't be 
necessary if the ultimate functional were known. 

(iii) Mixtures of 3 He and A He droplets. 

Mean-field like calculations with effective interac- 
tions are also used to describe those systems [13( 
and can be justified more fundamentally by the mul- 
ticomponent internal DFT formalism. 

The article is organized as follows. We first give the 
foundations of the multicomponent internal DFT formal- 



ism and make explicit the link with the previously de- 
velopped formalisms. Then, we apply the formalism to 
molecular systems where the nuclei are treated explic- 
itly. This application will show clearly the link between 
multicomponent internal DFT and "one kind of particle" 
standard and internal DFT. Finally, we mention some 
features of the application to atomic nuclei and mixtures 
of 3 He and 4 He droplets. 



II. MANY BODY FORMULATION. 
A. General formulation. 

For the sake of simplicity and to underline more clearly 
the physics, we assume: 

(i) Two different kinds of particles. 

This is sufficient to describe atomic nuclei, mixtures 
of 3 Hc and 4 He droplets and molecular systems with 
only one kind of nuclei. The generalization to more 
kinds of different particles, for the description of 
molecular systems with more than one kind of nu- 
clei, is easy. 

(ii) Two-body particle-particle interactions. 

The generalization to 3-body etc interactions is 
straightforward. 

(iii) Particles without spin. 

Generalization to spin polarized systems can be ob- 
tained by introduction of an arbitrary internal mag- 
netic field J2 i 'B int (ri — R), which models the ef- 
fect on internal properties of magnetic fields used 
in experiments, and adapting the standard deriva- 
tion pm. 

We thus consider: 
(i) iV^ 1 ) particles (1) of mass rrS 1 ^ and laboratory co- 



ordinates {r^ 



(i)J' 



(ii) iV^ 2 ) particles (2) of mass m^ and laboratory co- 



ordinates {r^ 
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and start from a general translationally invariant many- 
body Hamiltonian l 



H 



N m (1) - 
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1 In a Hamiltonian and wave function based description of an 
isolated self-bound system, the Hamiltonian should be explicitly 
translationally invariant to ensure Galilean invariance of the wave 
function. Translational invariance, which states that the observ- 
ables do not depend on the position of the cm., is a necessary 
but not sufficient condition for the more fundamental Galilean 
invariance, which ensures that scalar observables are the same in 
all incrtial frames. 
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composed of: 










(i) the usual kinetic energy 


terms, 
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(ii) translationally invariant two-body potentials u^ 1 ', 
u^ and u*- 12 -* which describe the particle-particle 
interactions, 

(iii) arbitrary translationally invariant " internal" poten- 
tials v\J t and v\ nt , which act on each species in the 
cm. frame and can model internal effects of poten- 
tials used in experiments (such as electric fields in 
the molecular case) , or can be safely dropped at the 
end in the case of an isolated self-bound system. 

The total cm. coordinate R, i.e. of particles (1) and 
(2), is defined as 
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where M^ l > = N^m^ is the total mass of the particles 
(I) and RO = -t^jj X)i=i r i i s the center of mass of the 
particles (I). The (7V (1) + N^ - 1) Jacobi coordinates 
{£ Q } are defined as (see appendix \X\ for some reminders) 



for ae[V,NM-l]: £ a = r^ - A £ rf } , (3) 
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NW 
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for a = iV (1) : 

for a £ [N {1) + 1; N {1} + N {2) - 1} 



&. 



(2) 
" a _jvW+l 



m(2) E «- J v« rf ) 
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(a - iVW)m( 2 ) + MO) 



The £ a are relative to the cm. of the other 1, . . . , a — 1 
particles and are independent of R. They are to be distin- 
guished from the (A^ 1 ) + A^ 2 )) "laboratory coordinates" 
rf\ and the (7V (1) + A^ 2 )) " c . m . f rame coordinates" 

(r\ — R) relative to the total c.m. . 

It is to be noted that, for multicomponent systems, 
those "standard" Jacobi coordinates cannot be associ- 
ated to one specific kind of particle, in contrary to the 
laboratory or c.m. frame coordinates. More precisely, 



even if the {£ a ;a G [1;A"W — I]} are the Jacobi co- 
ordinates associated to particles (1) only, the {£, a ;a € 
[iVW; ATW + A^ 2 ) — I]} mix the laboratory coordinates 
of particles (1) and (2). 

We mention that the "standard" Jacobi coordinates 
used in this paper do not constitute the unique set of 
coordinates which permit to separate the c.m. contribu- 
tion from the internal contribution in the Hamiltonian, 
as we will do thereafter. Another possible decomposition 
is (instead of ((3])) 

for 1 = 1,2 and a& [1;N {1) -I] : 

f(0_ r (0 _IVr (i) 

Sq — 1 a+l 2-^i i ' 

i=l 

and R< 12 ' = R< 2 > - R«, (4) 

in addition to the total c.m. coordinate R defined in 
Eq. ((5|). R( 12 ) is the relative motion of the c.m. of each 
species. The use of this set of coordinates permits one to 
introduce the Jacobi coordinates £a reduced to a given 
kind of particles (I), which is more symmetric. In this ar- 
ticle, we decided to use the set defined in Eq. (J3J because 
it permits us to obtain more directly the results of Ap- 
pendix [Bland £ )iV Al in the limit when one kind of particle 
has a much larger mass than the other kind, and does not 
change the final results nor the way to obtain them (the 
formulation in terms of the set ([Jl is straightforward). 

Since the {rf - rgj, {rf - rj m/i) } and {rf - R} 
can be rewritten as functions of the £ a , the interaction 
potentials u^\ u^ 2 \ u^ 12 \ and the internal potentials 

v\nt an d vlnt can b e rewritten as functions of the £ Q . We 
denote: 

N w 

5>«(r«-rf) -> V®&,>.,tNM + N»>-x) 

i>j 

Af(i) Af(2) 

i=i j=i 

N m 

5>&(rj'>-R) -+ V&(t 1 ,...,Z Nm+Nm _ 1 ). (5) 

i=l 

After having defined the conjugate momenta of R and £ Q 
(see appendix |A~|) . we can separate dTJ) into H = Hqm + 
H- lnt , where the c.m. Hamiltonian (M = N^nv- 1 ' + 
N^m^ is the total mass) 
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is a one-body operator acting in R space only, and the 
internal Hamiltonian 
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flint = E 
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+U '(£1, ...,^jv(i)+iv(2)-i) + U ( '{£,1,— Anw+nW-i) 

+ V int (£l>— i^VCO+ATW-l) + V int (£l) — > ^JVW+iVW-l) 

is a (A^ 1 ) + A"*- 2 ) — 1) body operator in the {£ Q } space. It 
contains the particle-particle interaction potentials and 
the internal potentials. r a is the conjugate momentum 
of £ a and /Xq, the corresponding reduced mass, defined as: 

a (1) 
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/or ae [7V (1) ;Ar (1) +7V (2) -1] : 

P + (a - iW) + l)m( 2 ) 
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Hence, the (A^ 1 ) + A"( 2 ))-body laboratory wave func- 



tion ip(r\ , . 



.(i) ._(2) 



, rU 2) ) can be separated 



into a wave function T that is an eigenstate of Hcm and 
depends on the cm. coordinate R only, and an "inter- 
nal" wave function ipi nt that is an eigenstate of H{ nt and 
depends on the remaining (A^ 1 ) + N 1 - 2 ' — 1) Jacobi coor- 
dinates £ Q : 



ip(r\ 
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r (l) , r (2) 
■ ' 1 AT(i) ' 1 1 ' ■ • 

r(R)Vm t (a,. 
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The r(R) describes the motion of the cm. of the isolated 
system in any chosen incrtial frame of reference (such as 
the laboratory). The ipint describes the internal prop- 
erties and is a function of the (N^> + N^ — 1) Jacobi 
coordinates. Of course it could also be written as a func- 
tion of the A^ 1 ) coordinates r\- and A^ 2 ' coordinates 
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rU 2) ), but one of 



them would be redundant [261 ] . 

Thus, cm. properties and internal properties arc fully 
decoupled and the total energy splits into E = Eqm + 
Ei nt . Since T is solution of the free Schrodinger equa- 
tion, r(R) should be an arbitrary stationary plane wave, 
i.e. infinitely spread and thus not normalizable. This 
leads to derealization of R and arbitrary cm. energy 
Ecm — ft 2 K 2 /(2A-f). This does not correspond to exper- 
imental situations, where the system is no longer isolated: 
interactions with other systems of the experimental ap- 
paratus localize the cm.. But the formal decoupling be- 
tween cm. motion and internal properties permits us to 
left the cm. motion to the choice of experimental condi- 
tions, internal properties being comparable to the exper- 
imental ones. If ipmt is normalizable, which is by defini- 
tion always the case for the ground state of a self-bound 
system, the internal energy can be written 



Ei„t[fpint] = (tpint\H in t\lpint)- 



(9) 



relations that will be useful for the subsequent consider- 
ations. Following Refs. [y, [l0|, [27|, we define the internal 
one-body densities for each species 



p W(r) = AT« ldr^...dr^LdA 2) 



l AT(l> 



■ dr 
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r (2) ^,2 



x8(K)\i> in t(r\ ; ,...,r Ar(1) ,x 1 , 

x5(r-(r«-R)). 

They are normalized to N"' . The laboratory densi- 
ties p"'(r) arc obtained by convolution of p in \ t with 
the cm. wave function (following [Hi, l27|): P^'( r ) = 
/dR|r(R)| 2 p«(r-R). 

Following the considerations of Refs. [y, [23], we define 
the local parts of the two-body internal density matrices 
for each species 

7&(ry) = N«HNV-l)JdrU...dTV ) drW...dr% 

x<5(R)|^ mt (rW,...,r« ); r( 2) ,...,r^ 2) )| 2 x 
x(r-(v^-Tl))8(r'-(r%-R)). (11) 

These have the required normalisation to N^-y (N® — 1). 
Following similar steps to those in Refs. [10, H3, we 
can show that the local part of the two-body labora- 
tory density matrices 7^(r, r') are obtained by convo- 
lution of ~f\ nt with the cm. wave function: 7®(r, r') = 

/dR|r(R)| 2 7 «(r-Rr'-R). 

Finally, we introduce the local part of the two-body 
internal "coupling" density matrix 

7&2W) (12) 

= N^N^ Jdr^...dv%dr^...dr% 

x8(R)\i }int (v^\...,v%-vf\...,v%\ 2) )\' 

x6(r - (r 2 (1) - R))<5(r' - (vf - R)), 

where r acts in the particle (1) space and r' acts in the 
particle (2) space. 

The definitions of p^ t (v), 7 £ ) t (r, r') and 7 ^ 2) (r,r') 
show clearly that they are defined in the cm. frame, 
i.e. that the positions r and r' are measured in the cm. 
frame (see the 8 relations in Eqs. (JTUJ), (fTTj) and ([12")) ). 
Compared to the traditional definitions, a 6 (R.) appears 
in the definition of the internal densities calculated 
with ip int written in terms of the {r\ } coordinates. 
Since one of them is redundant, the 8 (TV) represents 
the dependence of the redundant coordinate on the 
others 2 . Note that, following steps similar to those in 



B. Some useful definitions. 

Before we adress the generalization of internal DFT to 
many kinds of particles, we define some quantities and 



2 More generally, we should introduce a <5(R, — a), where a is an 
arbitrary translation vector, instead of the <5(R). This is linked 
to the translational invariance. In this paper we chose a = 0, 
which leads to perfectly equivalent results and permits to lighten 
the notations. 
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Ref. [6|, one can rewrite p\^ t , 7^ and ^l nt ' as functions 
of Jacobi coordinates {£ Q }. 

We establish a useful relation. For any func- 
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dinates expressible in terms of Jacobi coordinates [we 

denote F(£i, ■•■,Cjv< 1 )+A r < 2 >-i) an d -^ tne associated oper- 
ator] we have 
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We see that the "internal mean values" calculated with 
ipint expressed as a function of the [N^ 1 ' + iV^ 2 ) — 1) 
coordinates £ a , can also be calculated with ipi n t expressed 
as a function of the (A^ 1 ) +N^) coordinates r^ ' . Then, 
as above, a 5(R) that represents the dependence of the 
redundant coordinate on the others appears. 



with respect to the laboratory coordinates, becomes one 
body (and local) when expressed with the cm. frame 

coordinates (remember that the p\^ t are defined in the 
cm. frame, cf. JjIIBI) . 

The internal energy Ei nt [ipi n t] © can be rewritten 



III. GENERALIZATION OF INTERNAL DFT 

AND KOHN-SHAM SCHEME TO MANY KINDS 
OF PARTICLES. 

A. Hohenberg-Kohn theorem. 

We show that ipint can be written as a functional of the 
internal densities p\ n \ and p\ n \. The relation ([IS)) leads 
to (V£l is defined with ©) 



(iptntlVinll^nt) 
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As in its definition enter two arbitrary one-body poten- 
tials in the cm. frame of the form fdr vlJ t (r) p in \{r), and 
as the ground state ipint of Hi nt is obtained by minimiza- 
tion of Ei n t , we can directly apply the usual proof of the 
standard HK theorem [l|, [j] generalized to many kinds 
of particles [13, |2Q|- We thus can claim that, for a non- 
degenerate ground state ipint and given kinds of particles 
(1) and (2) (i.e. given interactions iS 1 ', vS 2 ^ and u^ 12 '), 
the ground state ipint can be expressed as a unique func- 
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also the corresponding internal energy of a self-bound 
system, i.e. E int [p\^ t , p\ n \]. 

As emphasized in jy, [30], the HK theorem is valid 
only for arbitrary one-body potentials that lead to bound 
ground states. As a direct consequence, the internal DFT 
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formalism is valid only for potentials v int and 
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where we used (|10[) to obtain the penultimate equality. 

,.(l) (r (l 

J mt K 1 1 



We see that the potential v ir }t( r i ' ~ R) that is iV^-body 



lead to bound internal ground states ipint- For pure self- 
bound systems, described by our formalism in the limit 

v int = v int = ^i ipint should by definition be a bound 
ground state, so that the previous conclusions still hold. 



Internal Kohn-Sham scheme. 



the interacting kinetic energy 
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Following similar steps as in (fT3|) . we rewrite simplicity) 
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which permits to recover an interpretation in terms of expression. 



particles (1) and (2) and makes it clear that the major 
difference with the standard kinetic energy comes from 
the cm. correlations, i.e. the <5(R) term in the previous 



Using relation (fTT)) and the notations introduced in 
£)IIB[ we can rewrite the internal energy as 



J 



P ^W n ( 2 )l 



dv^---drf 



N w (1)2 JV' 2 ' (2)2 



2m (i) ' ^ 2m< 2 

2—1 2—1 



\ j 'drdr' 7 ^(r,r') M (1) (r - r') + \ j "drdr> ^](r,r')u^(r - r') 
drdr , 7£f(r,r')u (ia >(r-r') 

^l(r)p«(r) + |dr^)(r)^(r). 

I 



r (2) 



>) 



(18) 



To recover the "internal" KS scheme, we assume, as in 
the traditional KS scheme, that there exist, in the cm. 
frame (the p\^ t being defined in the cm. frame, see 
^IIB[) . two non- interacting systems (i.e. two local single- 
particle potentials v s ): 

,2 



f-P- 



,c: 



V2m« 



(r)W 



1 = 



e (l) ^ l)l (r) 



1,2 , i = l...JV w ,(19) 

which reproduce exactly the densities /9^ nt of the inter- 
acting system: 



p&M 



jv(0 

E 



^Wl 5 



(20) 



In Eq. (fl9| we implicitly supposed that the particles 

are Fermions (i.e. that the {^n*} are orthonormal) . A 
KS scheme to describe Boson condensates can be defined 
similarly by choosing identical f^l for a given kind of 

particles (I). Uniqueness of the potential v s (r) for a 

given density p\^ t (r) is ensured by a direct application of 
traditional DFT formalism to each KS equation (|TT)|) . Of 
course, the question of the validity of the KS hypothesis, 
known as the non-interacting v-representability problem, 
remains, at least qualitatively, the same way that in tra- 
ditional DFT [l| (but can quantitatively be different). 

To use notation similar to the traditional no- 
tation, we add and substract from the inter- 
nal energy f| 18[) the internal Hartree energies 



Note that, even if only (N^ + N^ -I) coordinates are E H \P™t\ ~ 2 / '*" dr> Pint( T ) Pint( r> ) u r 'X 



sufficient to describe internal properties, we introduce in 
the KS scheme (A^ 1 ' + N^ 2 >) orbitals, i.e. the same num- 
ber that the number of particles. The reason is discussed 

in§|mcl 



the "coupling" Hartree term 

p (l2)r (1) (2) 1 ( drdr' n {1) (r\n {2) 

&H [PinUPintl — J arclr Pint\ r )Pint 

and the non-interacting kinetic 



iJV 



E i= i (v 



(0*i 



p 2 i..(0<\ 



( r ') u (12)( r _ r /) 

energy terms 
This permits us to rewrite 



the internal energy (|18|) as 



(1) (2)l 
irit\Pinti Pint! 



Ei n t [p. 



(i)i 

™* l 2m( 1 ) 



(i)i 



lyl^int) 



(2)r (2)-, 



W&l + ^Mni] 
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2m( 2 ) 

P (12)r (1) „(2)- 
"ff [Pint' Pint. 



I (2)i\ 



+£ 



,( 2 )l 



XC\Pinti Pint! ' ^XCVP 



<2) 



■/■ 



,(D 



.Jl 



+ /drt/£J(r)p&(r)+ dr^(r)p^(r). 



,(2), 



, p( 12 )r /) (l) „(2)i 
J+^C [Pint 'Pint J 



.(2), 



(21) 



The internal exchange-correlation energy for the iden- 
tical particles of kind (I) is defined as 



P W r„(i) n (2)i 

^XClPinf Pint! 



2 /rfrdr' 



Aft" 



7 «(r,r')-p^(r)p«(r') 



uWfr-r'l 



,(2) 
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( 2 )) 2 m(0^ m ^ ri 



(i) 



(22) 



,(2) 
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)~( ( /'intl^Z7nl¥'ilt) 



,(0< 
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Note that, since the KS assumption implies (p\ 7 il[pj lt 



(i) 

in I 



a, 



and p^ t . 



-5^(7 can we ^ ^ e wr hten as a functional of p. 



We see that E xc contains the exchange-correlations 
that come from the interaction u"' (first line of (|22|) ). but 
also the correlations contained in the interacting kinetic 
energy (second line of (|22|)). Concerning these correla- 
tions, it is clear that they come, on the one hand, from 
the correlations neglected in the traditional independent- 
particle framework, but also from the cm. correlations 
(the 5 (TV) term in the previous expression). The kinetic 
energy term is the only one that explicitly contains those 
correlations because they directly affect the motions of 
the particles in the cm. frame. Indeed, in this frame, 
if one particle moves in one direction, the other particles 
will tend to move in the opposite direction. As quantum 
particles are described by wave functions, they always 
have associated zero-point motions. Those zero-point 
motions are coupled by the cm. correlations, which is 
specific to quantum systems; the corresponding quantum 
part of the cm. correlations is included in E xc 3 . 



3 We take the opportunity to underline a key difference between 
the classical and quantum separation of the cm. motion. In 
classical mechanics, as particles are pointlike, this separation is 
done simply by a coordinate change. It can be done using Jacobi 
coordinates (the redundant coordinate is then treated explicitly) 
or the cm. frame coordinates (which lead to a similar final re- 
sult because the redundant coordinate is implicitly taken into ac- 
count through the translational symmetry of the Hamiltonian). 
In quantum mechanics, particles are described by wave functions. 
Then, one has to separate the non normalizable part of the lab- 
oratory wave function, which can be done using Jacobi coordi- 
nates. Another consequence is that, in the cm. frame, all points 

of space that satisfy m' 1 ' Y.?=i r f 3 + ™ (2) Efci' r ?' = ° ( and 
only those points) are allowed. This couples the zero-point mo- 
tions and produces a purely quantum contribution to the cm. 
correlations. 



The inclusion of the cm. correlations in the exchange- 
correlation functional is the major difference with tradi- 
tional DFT 0, Q and previously developped multicom- 
ponent DFT formalisms @, Il6l - t20| , and opens the way to 
the search for a local cm. correlations potential, which 
would a priori be computationally much less costly than 
the projection techniques used, for instance, in nuclear 
physics [l2, I21H23J . Note also that the inclusion of the 
cm. correlations explains why E xc is a functional of 



,U) 



.,(2) 



(I) 



both p\ n ' t and p\ nt , and not only of p\^ t as one might 
have expected: those correlations necessarilly couple the 
particles (1) and (2), thus their densities (this will be 
underlined from another point of view in £)III Dp . 

We emphasize the fact that a part of the total cm. 
correlations is contained in E xc and another part in 

E xc . Indeed, it is the interacting kinetic energy term 
as a whole that contains explicitly the cm. correlations, 
and this term is split in both E xc and E xc . The decom- 
position we choose for the energy functional is of course 
not unique. But even if it does not include all the cm. 
correlations in one specific functional, its advantage is 
that the meaning of the obtained functionals is clear and 
that it permits us to recover the functional forms of " one 
kind of particle" traditional and internal DFT in the lim- 
its that will be discussed in i jIVl 

The internal coupling correlation energy between par- 
ticles of kinds (1) and (2) is defined as 



j(i2) r „(i) „(2)i 
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En ' I Pint 1 Pint} 



(23) 



„(!2) 



drdr' 7 ^(r,r')-p^(r)p^(r') 



(i) 



.(2) 



U ( 12 )(r-r' 



There is of course no exchange energy in (f2"3"|) because the 
particles (1) and (2) are not identical. As this term does 
not contain explicitly the cm. correlations, the spirit 
of the already developped approximations found in [8|, 
llq . |2C§ remains suitable within the formalism presented 
here, at least for molecular systems, especially that of [16[ 



which is directly applicable to a "coupling" correlation 
energy written in the form (f2U)l . 

The remaining task is to minimize the internal energy 
(f2"Tj) so as obtain the equations of motion (which define 

] with respect to <p\„l* ' , and 



Pint ) ■ Varying E mt [p^ t , p[l\ 

imposing (ortho-)normality of the {f\„l} 



leads to "internal" Kohn-Sham equations (fT9|) for the 
Wilt) with 






(r) 



[E ln t [p 



(i) (2)1 

inf I int\ 



2 iV (m) 
?n— 1 z— 1 
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= [dr'pW(r')u«Hr - r') + !7&V&>p£l](r) + «&(') 



drV&V)^ - r') + E#*"V&P&](*) + ET^r-W ,(2) 



Pint. Pint! ( f ) 



/or m ^ L (24) 
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We keep v znt for generality but recall that it can be set to 
zero in the case of an isolated self-bound system. In addi- 
tion to the classical parts of the interaction (the Hartree 
potentials for the (I) — (I) and (1) — (2) interactions, and 
the optional external potentials) appear purely quantum 
potentials, which are defined as follows. 

The exchange-correlation potential for the particles of 
kind (/) is 



jcj?(0 r (!) ( 2 )i 
r/C/Or-fi) n (*h( r \- 5E xc\PinvPint\ 

U XC IPinti PintlV 1 ) — _. (ft . . 



,(!) 



J 2 ) 



(25) 



(') 



It is functional of p int and p int , and not only of p\ nt , 
for the reasons discussed above and linked to the cm. 
correlations. This potential contains the "standard" 
exchange-correlation for the (I) — (I) interaction and the 

part of the cm. correlations that is contained in E xc . 

The other part of the cm. correlations is included in 
^xc (because of the decomposition we choose). Thus 
there also appear complementary pure cm. correlations 
potentials 



U. 



im/l) [p^P^\(r) = 



(m) 



*W6\p 



J int' Pint! 



Sp&(r) 



for m 7^ I. 



Finally, the coupling correlation potential is 



ug 2/,) ip, 



inti Pint \ \ r ) 



dtj C [PinvPintl 



sp$ t <?) 



(26) 



All those potentials are local as expected, which is 
the computational power of DFT. One can note that 
Eqs. (fiT)]) and (j2"4")l are symmetrical under the exchange 
(1) o (2). 

Note that the first line of Eq. (f2"4"]l contains the terms 
that describe the interaction between the identical par- 
ticles of kind (I), which are close to those found in the 



"one kind of particle" KS scheme [5(- The second line 
of Eq. ((24)) contains terms that are specific to the KS 
equations for two kinds of particles. They describe the 
interaction of the particles of kind (I) with the particles 
of kind (m ^ /). 



C. Why do we introduce a number of orbitals 
equal to the number of particles in the KS scheme? 

The internal DFT formalism raises the following ques- 
tion (we consider only one kind of fermion in this section 
to simplify the discussion): since, for a self-bound sys- 
tem composed of N particles, only (N — 1) coordinates 
are sufficient to describe internal properties, why do we 
introduce N orbitals in the internal KS scheme? We men- 
tion that the question of the number of orbitals to intro- 
duce in the KS scheme also occurs in traditional DFT. 
Indeed, this number is not imposed by the theory; it has 
been fixed "by hand" to N orthonormal orbitals (intro- 
ducing an auxiliary system of N non-interacting parti- 
cles). From now, nothing would fundamentally forbid 
us from introducing another number of orbitals (eventu- 
ally non-orthogonal with effective masses). This would 
of course change the form of the KS energy functional, 
which will "adapt" to the change of the number of or- 
bitals. More precisely, since to a choosen number M 
of orbitals corresponds a specific definition of the non- 
interacting "kinetic energy", it would change the part of 
Exc[p\ that corresponds to the difference between the 
interacting kinetic energy and the non- interacting "ki- 
netic energy" (i.e. the second line of ([22|) ). that we note 
EAkin[p\- For instance, in the case where only one KS 



orbital would be introduced, E 



(M=l)r 



^Afcm [PI WOuld bc draS " 

'[p], but it is easy to show 
that the corresponding single KS equation would lead 



tically different from E ( Akin '' 



to the exact interacting ground-state density, i.e. to the 
same result that the one obtained by direct variation of 
the HK functional by p, so that the non-interacting v- 
representability is perfectly achieved (since the KS non- 
interacting " kinetic energy" can then trivially be written 
as an explicit functional of p ((pi = ^fp) and thus its vari- 
ation is defined for all densities for which J~p is differen- 
tiable, including the interacting v-representable ones). 

Purely formally speaking, the question of the optimum 
number of orbitals to introduce in the KS scheme can 
be reformulated as follows: is the non-interacting v- 
representability better achieved with a certain number 
of orbitals M (> 1)? This question is still open in tradi- 
tional (and internal) DFT. 

But practically speaking, many advantages justify the 
introduction of N orbitals: 

(i) Even if only (N — 1) coordinates are sufficient to 
describe the internal properties of a self-bound sys- 
tem, they still describe a system of N particles. We 
thus have to introduce N orthonormal orbitals if 
we want them to be interpreted (to first order only) 
as single-particle orbitals and obtain a scheme com- 
parable to actual mean-field like calculations with 
effective interactions. In particular, the highest KS 
eigenvalue e^r is then the negative of the ionization 
(or the separation) energy with the exact functional 

HE HI. 

(ii) The anti-symmetrization is explicit in terms of the 
N particles coordinates only. Introduction of N or- 
bitals permits us to explicitly take into account the 
anti-symmetrization in the definition of p, while the 
anti-symmetrization would be implicit if another 
number of orbitals was introduced and thus more 
difficult to obtain because it would have to be more 
precisely taken into account in parametrizations of 
Exc[p\- 

(iii) It leads straightforwardly to the classical (point- 
like) limit, see appendix iDl 

(iv) It a priori permits one to obtain, in the general case, 



the non-interacting kinetic energy that is the clos- 
est to the interacting kinetic energy (i.e. the small- 



J Akin l 



est E\j„\fi\ 



Indeed, E^^ [p] contains only the 
correlations energy part of the interacting kinetic 
energy, which is generally a correction. It is essen- 
tially for quite small self-bound systems and within 



(M=N) 
Akin 



[P\ 



can 



the internal DFT formalism that E 
be more than a correction, because it contains the 
cm. correlations, which can be large for those sys- 
tems. Then, one track is to conserve N orbitals and 
to develop precise parametrizations for E^ k ~ n [p] 
(which is not achieved by actual functionals); this 
would permit us to preserve the advantages listed 
above. But the previous discussion also let the door 
open to another possible track: study whether the 
introduction of a number of orbitals M =^ N (for 
instance M — N — 1 with effective masses) would 
permit us to obtain a smaller E Ak } n [p] , so that dif- 
ficulties in parametrizing E^ k - n [p] would affect the 



result as little as possible, 
scope of this paper. 



This goes beyond the 



D. Link with previous multicomponent DFT 
formalisms. 



The previous multicomponent DFT formalisms found 
in Rcfs. [16- 19] have been developpcd without separation 
of the cm. motion, as already underlined in Ref. (18J . 
The definitions of the energy functionals were based on 
the laboratory wave function ip and densities p"> . To 
understand the key differences with those previous for- 
malisms, we switch to the Levy-Lieb constrained search 
formulation [30, |31j , within which they have been devel- 
oppcd. The corresponding energy can be decomposed in 
a form similar to Eq. ([2~T]) . where only the definitions of 



Ey n and E]-, change (i.e. the fourth line of Eq. (f2"Tjl ) 
[Sill El: 



ExClPint' Pint! ~~ ^ 

E%[ P V] 

4 12 V 1} ,p (2) ] 



j/iCO-^TO 



(27) 



(^OlfW + tfCO^CO). min ^f D \fV\i,f D )-E$\p% 



min (MT^+T^ + U^ + U^+U^\^)- min (V> (1) |T« + U^ k/>«) 

l/>-S-p (1) ,p (2) ^(l)^pd) 

- min (^ 2 )|f( 2 ) + LX 2 )|^ 2) )-^ 12) [p (1) ,p (2) ]. 

lfj( 2 ) _).p( 2 ) 



r 



The if; are the interacting (iV^ + iV' 2 ))-particles states, (Slater Determinants), the ip^ are the interacting NW- 
the ipg D are the non-interacting N® -particles states 
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(!) p(l)2 



is the kinetic en- 



particles states and T"> = ~J2 i=1 — — — 
ergy operator for particles of kind (I). The interest of 
the introduction of the tp (and of the constrained search 
formulation) is that it permits us to define E xc as a func- 
tional of p( 1 ' only (and not of both pW and p^)- The 
constrained search formulation is by definition restricted 
to densities and thus states that are normalizable (the 
variational principle being defined only for systems for 
which the wave function is normalizable [25, 30]). How- 
ever, as the functionals are formulated in terms of the 
laboratory density (i.e. the cm. motion is not sepa- 
rated) , the equations obtained after minimization of the 



energy functional should lead to full derealization of the 
density. Moreover, as the theory is not formulated in 
terms of internal properties, the cm. correlations do not 
appear explicitly in the functional. 

Internal DFT overcomes those problems because it is 
based on the internal wave function, which is by defini- 
tion always normalizable for self-bound systems, and is 
the one which defines the observables of experimental in- 
terest [6[ . For better comparison with the previous multi- 
component DFT formalisms, we reformulate the energies 
E x ' c and E^, defined in i]IIIBI within the constrained 
search formulation: 



J 



P W r„U) „(2h 



(12) r (1) (2), 



E'c'lP 



mm 

/ , (1) (2) 

»/W->- P}n t , Pint 



I , (1) (2) 



ftMrW + U^int) ~ min (4'L|f W|4^) - E$\p® 

Vsn^Pint 

(^„t|f« +fW +#W +fj^ +W 1% Mint) 



(28) 



, . (1) (2) 



(^nt|T (1) + L>W|^„t) 



/ . (1) (2) 



(An t \f^+U^ lnt )^E^[p^ t ,p^ t }. 



The tpint are the interacting internal states which can be 
expressed with (N^ + N^ — 1) Jacobi coordinates 4 . 
The ip SD are the non-interacting iV^-particles states. 

Note that the choice we made in (|28p for the ip used 
in the minimization of (V>|T^ + tA"|V>) is not unique 
(because there is a cancelation between E xc and Eq ) . 
One could, for instance, use interacting states ip, t ex- 
pressible with (TV") — 1) Jacobi coordinates which lead 
to p i7 l t , instead of the (N^+N^ — l) Jacobi coordinates 

state leading to p int and p int 5 . Even if formally correct, 
this is not the choice used in (|28|) mainly because it would 
lead to interpretation difficulties. Indeed, it amounts to 
considering that the particles of kinds (1) and (2) have 
independent cm. motions 6 , while the cm. motion is 
common for the system as a whole. Moreover, this would 
bring no gain because it would move the explicit inclusion 



of total cm. correlations to E, 



(12) 

c 



so that previously de- 
tcrm dGl [13 wou ld 
not be fully usable, contrary to the formulation proposed 



lat previoui 
velopped approximations for this term [8j, lla . |20| would 



4 Note that when the integrals involving ipi n t are written in the 
(TV' 1 ' +JV( 2 ') particles coordinates representation, following sim- 
ilar steps as in 11311 , a <5(R) appears explicitly. 

5 Since the cm. motion is "substracted" in our formalism, the 
most pertinent interacting states to use are of course internal 
ones, i.e. expressible in terms of Jacobi coordinates. 

6 Indeed, following similar steps as in 1131 1. the integral {w^^T "' + 
U"'\i/>.* t ") can be written in terms of the JV"' coordinates of 
particles (I). Then, a <5(R.'") appears explicitly. 



in this paper. Another advantage of the decomposition 
we choose in this paper is that it permits a clear link with 
"one kind of particle" standard and internal DFT, which 
will be discussed in mV\ 

We now discuss the link with the slightly different work 
found in |8|, [2(| ■ First of all, this work is by its very na- 
ture essentially suited for the description of molecular 
systems, since the electronic coordinates are transformed 
to the body-fixed frame independently from the nuclear 
coordinates (this "decoupling" is possible in molecular 
systems due to the large difference of masses between 
electrons and nuclei, see appendix [B] independent trans- 
formations arc no longer possible for a self-bound system 
where there is no large difference of masses between the 
different kinds of particles constituting it). Second, this 
work treats the breaking of rotational symmetry, which 
goes further than the other multicomponent DFT for- 
malisms (we don't treat rotational symmetry in this ar- 
ticle, but we are working on the problem). Third, the 
nuclear iV^-body density is used as a basic variable in- 
stead of the nuclear one-body density, for easier descrip- 
tion of collective phenomena (but it is of course less easy 
to handle numerically for quite large systems). 

The difficulty discussed above remains: even if the 
electronic coordinates are transformed to the body-fixed 
frame, since the cm. wave function is not separated from 
the beginning, all the densities are defined through a fully 
delocalized wave function (for a ground state). But nu- 
merically speaking this does not cause a real problem 
because: 

(i) The KS step breaks by essence the translational 
symmetry, so that it forces the densities to become 
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localized. 

(ii) In Refs. [||,[2(j] the KS scheme is formulated in terms 
of the 7V(™)-body nuclear density and, in the pre- 
sented practical calculations, the nuclear cm. has 
been separated at a second step, leading to the in- 
clusion of the cm. correlations in the corresponding 
N^-body auxiliary nuclear wave function. 

In the formalism presented in the present article, the 
cm. motion has been separated from the beginning, lead- 
ing to a proper formulation of the HK theorem in terms 
of the internal densities and the explicit inclusion of the 
cm. correlations in the functional. This formalism is 
suited to study all self-bound systems. 

We now study thoroughly the application to molecu- 
lar systems, which permits us to make explicit the link 
between this formalism and traditional DFT, and then 
detail the application to atomic nuclei and mixtures of 
3 Hc and 4 Hc droplets. 



IV. APPLICATION TO MOLECULAR 
SYSTEMS AND LIMIT OF TRADITIONAL DFT. 



the nuclei will be described by (N^ — 1) Jacobi coor- 
dinates and will carry all the cm. correlations. As a 
consequence, the electrons are not concerned by the re- 
dundant coordinates problem (i.e. the cm. correlations) 
and remain described by N^ e ' coordinates in the frame 
attached to the cm. of the nuclei. Sec appendix [B] for 
better understanding of this mechanism from the point 
of view of Jacobi coordinates and appendix [C] from the 
point of view of the functionals. 

As a direct consequence, the electronic exchange- 
correlation functional E xc no longer contains explicitly 
the cm. correlations and its form becomes comparable 
to that of the exchange-correlation functional of "one 
kind of particle" traditional DFT [l[, see Eq. JC3]). All 
the cm. correlations are explicitly included in the nu- 
clear exchange-correlation functional Ej?q only and its 
form becomes comparable to the form of the exchange- 
correlation functional of "one kind of particle" internal 
DFT @, see Eq. JC4]). Since it is only the cm. cou- 
pling between the particles (I) and (to 7^ I) that produces 

a p^ dependence in E^' c (see discussion following 
Eq. ([22])). we have: 






for m 7^ I 



Electrons (we note I = e) arc Fcrmions, so that ex- 
change and correlations arc included in the functional 
E^q. Nuclei (we note I = n) are approximated as point- 
like and can be Fermions or Bosons (in the latter case, 

(n) 

the functional E y x £ contains only correlations and all the The internal energy (J3TJ) thus becomes 
ifint are identical). The internal potentials uj nt can then 
describe internal effects of a polarization potential (i.e. a 
voltage applied to the system), for instance in the station- 
ary case, or internal effects of lasers used in molecular ir- 
radiation experiments, for instance in the time-dependent 
case 1331. 



E (l) r (0 



^int[P mt , Pintl 



N< 



E(* 






N<- 



E&> 



(e)h_P 



*>&': 



Simplification of the formalism due to the large 
difference of masses. 



In molecular systems, the nuclei are much heavier than 
the electrons, i.e. mS n > » m^ . It thus should be a very 
good approximation to consider R = W~ n > and apply Ja- 
cobi coordinates to the nuclear coordinates only, so that 



+EP [,£)] + E (°) [pW] + E (ne) [p (n) p « ] 

+EP c ip%) + E¥ c [p%]+E£% ( £,p% t } 

+ J drv£> (r)pt| (r) + J drv% (r)p£> (r) , ( 29 ) 



and the KS potentials (jIMf for the electrons and nuclei 
become 



J 



< e) (r) = [dr> P % t (r>)u^(v r') + U^\p^](r) +„«(,) + f dr' p^(r')u^ (v - r') + ^"V&pSK*). 



.» 



(r) = / dr' P ^ t (r')u^(r r') + U { ^\p^ t ]{v) + «W(r) + / dr' p£ t (r')u^(r r') + U { ^ ] [p% t , p^ t ]{v) . 



(30) 



Since E^' c [pf^ t ] is comparable to the exchange- 



correlation functional of traditional DFT, standard ap- 
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proximations remain pertinent, for example the widely 
used Local Density Approximation (LDA) 7 , see e.g. 
[33| . or its extension to the Generalized Gradient Ap- 
proximation (GGA) [34| . 

It remains an open problem to find satisfying numer- 
ically manageable approximations for E x ^[p^.], which 
among others would contains the cm. correlations en- 
ergy. It would permit us to study molecular systems 
where the quantum nature of the nuclei plays an im- 
portant role(for example, small molecules 14 and solid 
hydrogen [HI])- The LDA could not be used a priori 
to approximate it because the nuclei are generally lo- 
calized. It would probably be more pertinent to start 
from an approximation which exploits localization (and 
not derealization as LDA docs). A crude approximation 
would be to do the maximum localization approximation 
and replace, in the obtained functional, the correspond- 
ing classical density by the true density of the system. 
Appendix [D] details how to make properly the classical 
(pointlike) approximation for self-bound systems, which 
is not completely trivial because of the translational sym- 
metry. Then, the "kinetic energy correlations" part of 
Exc ( wmcn contains explicitly the cm. correlations) 
desappears and 

7 ^(r,r')-^(i W) 

see appendix |DJ We thus ob- 



use the semi-classical approximation |43| instead of the 
crude classical (point-like) approximation. 

Concerning E^' [pi^ t , p\ nt ] i we recall that the spirit of 
the already developped approximations found in Refs. [H, 
llq . [20|| remain suitable within the formalism presented 
here, especially that of Ref. [Tg|, but futher developments 
are certainly desirable. 

Finally, one can note that the KS equations (fTTJ)) and 
(|30[) still remain symmetric under the exchange (e) O 

(n), but the forms of E xc an( i E xc (thus of U X q" and 

(e/ ) . 

U X q ) are different: one contains explicitly the cm. cor- 
relations whereas the other one does not. 

We now make further approximations that permit us 
to recover traditional DFT and thus to understand better 
the approximations inherent to this formalism and how 
to eventually improve its results. 



Further approximations that enable us to 
recover traditional DFT. 



2^,i=\ Pint \ L )Pint \ l ) i 



where pQ? 



(n)i 1 9 
\<Pint 



tain E 



(n) . _^7V<"> p (n)r (n)i 
XC ~ 2-*ri=l ^H [Pint 



the self-interaction correction to the {n)-{n) Hartree en- 
ergy, as noticed in Refs. [18|,[35| for the standard case, and 



leads to the potential U X q" 



S n )i 



1. Electron-nuclear correlation energy neglected. 

We suppose that, in the total internal energy (1291) . 
the electron-nuclear correlation energy E x ^ is negligi- 
ble compared to the other energies. This in particular 
implies that \Eq | is much smaller than |SjJ^|. It is 
likely that this is a fair approximation for a certain class 
of molecular systems, such as: 

/ dr 'ptnt (r')u (n) (r-r') (j) £ W should at least conta i n the (n)-(n) Hartree self 



which corresponds to 



when applied to the state ip int . This approximation gives 
a single densities dependence and not an explicit total 
density dependence (this dependence is nevertheless im- 
plicit as soon as a local KS potential is imposed, which 
is the aim of the Optimized Effective Potential method 
[3a - l38j . because it implies <p\™t [p|™t])- This necessitates 
extra (numerically more costly) caution to ensure the or- 
thonormality of the {<^)™j }, see for example Refs. |39r - 

m. 

Even if this "maximally localized" approximation 
should be used a minima in the absence of other approx- 
imations, its range of validity is not clear. Indeed, the 
width of the nuclear single densities is not systematically 
very small [iME^l- Moreover, this approximation does 
not permit us to evaluate the cm. correlations contained 
in E xc (because they disapear in the classical limit, sec 
appendix [D] and footnote [3]). Thus, the question of a 

(n) 

suitable approximation for E xc as a functional of the 
one-body total nuclear density, including the cm. corre- 
lations, remains open. A possible approach would be to 



interaction (as discussed above), 

(ii) electrons and nuclei are different kinds of particles, 

(iii) nuclei are generally localized (thus not too far from 
classical particles). 

However, the range of validity of this approximation is 
not presently perfectly clear. Indeed, there exist molecu- 
lar systems for which the electrons-nuclei correlations are 
important for the description of some physical properties 
[l6|. Then, one has to remain at the previous step. 

Within this approximation, the total internal energy 
becomes 



Eint\p\ 



(«) (e) 

Hnt [Pint ' Pint 



(n)i { _P | 

^ x ^ mt ' 2m(") ' 



N M 2 

(n)i\ V~^ / ( e )*| P 



l ) + E^- 



; = 1 



2m( e ) 



\<P 



(e)t 



, p(")[„( n )i , p( e )r J e )i I p(" e )rJ n ) „( e )i 

+E H [Pintl + h H [Pintl + h H [Pint > Pintl 



+I%h\p&] + E% ; \p&] 



+ / dTv£l(v)ptl(r) + / drv&(r)p£ t (T), (31) 



7 I.e. maximum derealization approximation (Fermi gas) and 
replacement, in the obtained functional, of the Fermi gas density 
by the true density of the system. 



and the KS potentials (|30| become 



4 e) (r) = ldv'p^ t {v')u^{v-v') + U^\pt n \}{v) 
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( n ) / \ 
•J S ( r ) 



+*&(*) + j dv'p^{v')u^\v - r'), (32) 



We see that, for the electrons, we have recovered the 

traditional " one kind of particule" KS potential (as U x q 
does not contain the cm. correlations), with, in standard 
notation [lHa] 

v ex \r) = «W(r) + |dr'p$(r> (en) (r - r'). (33) 

The potential v ext contains the arbitrary one-body inter- 
nal potential v^ t and the Hartree part of the electron- <J( p _ p\ ncl) ), so that p^(r) -> £ 4 =7 S ( r - r*™^) and 
nuclear interaction (which is a functional of the nuclear 



2. Classical (point-like) nuclei. 

A quantum treatment of the nuclei may be necessary 
for an accurate description of small molecules [14J and 
certain molecular systems of intermediate size. For larger 
molecular systems, a classical treatment of the nuclei is 
generally sufficient because the nuclei become more and 
more localized in {r} and {p} spaces. A further ap- 
proximation is to suppose that the nuclei are perfectly 
localized in both spaces. The appendix [D] shows how 
to make this approximation properly in the case of self- 
bound systems. We note {r\ nc ,p\ nc } the positions 
and momenta of the classical nuclei in the cm. frame, 



V ^JV ( " ) (ncl) 

i.e. satisfying 2^ i=1 v\ 



We pose |^(r)| 2 

(ncl). 



and X)i=i P* = 0. 



<5(r 



{ncl), 



and l^&P)| 2 

(ncl) -y 



» 



,(e) 



PintiP) -^ Sili <Kp - Pi C )• The total internal en- 



one-body density p\ J ) . v ex is thus functional of v int and ergy ^ becomes ( rem ember that only a Hartree self- 



Pj-™j, which is allowed by traditional DFT (traditional 
DFT is valid for any v ext that is not a functional of the 
electronic wave function, thus density, and that leads to 
a bound state). The potential v ext , which is internal for 
(self-bound) molecular systems, becomes external for the 
pure electronic problem. This permits us to understand 
why traditional DFT is particularly well suited for the 
description of the electrons (only) in a molecular system, 
in the frame attached to the cm. of the nuclei. 

Note that, contrary to what is sometimes thought, 
traditional DFT does not, strictly speaking, involve the 
clamped nuclei approximation (and thus not the Born- 
Oppcnhcimcr approximation either). Indeed, traditional 

KS equations (from which follow the <p\nt) do n0 ^ nec " 
cssarily depend paramctrically on the nuclear positions. 
One is allowed to use an external potential of the form 
(f3"3"|). so that the nuclei are described with a spatial 
width (i.e. quantum- mechanically) , and to specify an- 
other equation that defines p\2 t ■ It is when this last equa- 
tion is not specified (for instance if we simply impose that 
each nuclei is represented by a gaussian of fixed width) or 
when the nuclei are treated classically (see next section) 
that the total molecular energy should be calculated for 
various nuclear configurations (various mean positions of 
the gaussians in the first case) to deduce the configura- 
tion of minimum energy. In the most general case, the 
considerations of this section permit one to understand 
that, from the point of view of the KS potential, tra- 
ditional DFT fundamentally requires that the electron- 
nuclear correlation energy E^c [Pint > Pintl * s neglected, 
but not necessarily the clamped nuclei approximation. 

Concerning the nuclear potential Vg , note that it is 
similar to the "one kind of particule" internal KS po- 
tential (as cm. correlations are included in E^c), with 
(using the notations of jg|) 



interaction correction remains in E^ c at the classical 
pointlike limit, see ^IV A[ and that in the stationary case 
each p\ nc ' should be null) 

^nt[P in t,\ r i J I — 



JV (e) 2 

E, (e)i,_P_ 
Wint I 2TO (e) 



Vint ) + ^H [Pintl + ^XCVPintA 



+ J drv£ t (r)p% t (r) + £ J dvp^ t {v)u^\v 






(n)l Y (ncl) _ (nd)x 



N 1 - 



7 j V int 



in) / (ncl) > 



(r 



(ncl) \ 

r , ) 



(34) 



and the electronic KS potential (|32|) becomes 

4 e) M = Jdv'p^ t {v')u^{v-v') + U^\p ( ti\(j) 



*&M 



J2 « (en) 

i=\ 



(r 



(ncl) 



(35) 



which now depends parametrically on the nuclear posi- 
tions {r-" c '}. Thus, the {<p,-„/} and pf^ t also depend 
paramctrically on the {r s }. We recover traditional 
" one kind of particle" DFT equation [lJ-J5|] with a classi- 
cal nuclear background: 



v ext (r)=vt ) t (r) + 



J2 u^ 



(r 



(ncl) ^ 



l (r) 



&*) 



dr'p^(r>( e ">(r 



The nuclear KS equations become obsolete in the sta- 
tionary case. Indeed, classical pointlike nuclei have no 
zero point motions, so that they fall to the bottom of the 
potential well, and the ground state of the whole molecule 
can be found by minimization of the total internal energy 
([Ml) for various nuclear configurations {r\ nc }. 



3. A method to improve the results of traditional DFT. 

The previous considerations permit us to set up a 
method to improve the results of traditional DFT in the 
molecular case: 

(i) Start from the classical nuclei approximation de- 
scribed in £|IVB 21 with electrons described by " one 
kind of particle" DFT and parametrized by the nu- 
clear positions. This is what is done in most of 
practical calculations. This is essentially suited for 
molecular systems whose nuclei are very localized, 
i.e. quite large molecular systems. 

(ii) If this is not enough, for instance in the case of rela- 
tively small molecules where quantum effects associ- 
ated with the nuclei play a role, skip to the approx- 
imation described in ^IVB 11 The nuclei are then 
treated quantum mechanically but the electron- 
nuclear correlation energy is neglected. There re- 
main the (n)-(n) Hartree interaction with at least a 
self-interaction correction, until other more satisfac- 
tory nuclear exchange-correlation functionals (in- 
cluding the cm. correlations) arc available. The 
electrons still satisfy "one kind of particle" tradi- 
tional KS equations, but the corresponding orbitals 
are no longer parametrized by the nuclear positions. 

(iii) If this is still not enough, for instance in cases where 
the electron-nuclei correlations play an important 
role, add the electron-nuclear correlation energy as 
described in JjIVAl using for instance the approxi- 
mation proposed in [16l | (but further devclopemcnts 
appear desirable). The obtained equations go be- 
yond "one kind of particle" traditional and internal 
KS equations, because of the coupling correlation 
term 4"%^ » 
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ternal densities of the protons and neutrons 8 . It gives 
a more fundamental justification than the Hartree-Fock 
(HF) framework to the stationary nuclear mean-field like 
calculations. Indeed, HF does not contain quantum cor- 
relations, nor does it treat correctly the redundant co- 
ordinate problem, which introduces a spurious coupling 
between internal properties and cm. motion [2ll. |44||. A 
way to overcome this problem in the stationary case is 
to perform projected HF (projection before variation on 
cm. momentum), which permits to restore Galilean in- 
variance, but at the price of abandoning the independent- 
particle description [12|, l2l| - |23| and at a large numerical 
cost. The internal DFT formalism demonstrates that the 
cm. correlations can be taken into account through a lo- 
cal potential, a priori numerically much less costly than 
the projection techniques. 

We underline that, for nuclear systems, the v\J t are 
generally zero. It is the dependence on the initial state 
that allows, in the time-dependent case, to describe for 
instance the collision of two nuclei in the frame attached 
to the total cm. of the nuclei [32| . 

We mention that the point of view often adopted in 
nuclear physics is that the protons and neutrons are the 
same kinds of particles, leading to isospin considerations 
[44J . From this point of view, it is rather the inclusion of 
the isospin in "one kind of particle" internal DFT that 
would give KS equations comparable to nuclear mean- 
field like calculations. 

The application of multicomponent internal DFT to 
mixtures of 3 He and 4 Hc droplets 9 also permits us to give 
a fundamental justification to the mean-field like calcula- 
tions done to describe those systems [13]. Since the cm. 
correlations are, to our knowledge, still not treated in 
existing calculations of Helium droplets, this work opens 
the way to their inclusion. 



VI. CONCLUSION. 



V. ON THE APPLICATION TO ATOMIC 

NUCLEI AND MIXTURES OF 3 HE AND 4 HE 

DROPLETS. 



When the masses of each species constituting the self- 
bound system are not very different (for instance in the 
case of protons and neutrons in an atomic nucleus or 
mixtures of 3 He and 4 He droplets), one cannot do the ap- 
proximations presented in the previous section, because 
the cm. correlations couple all the particles, and the 
complete formalism presented in J jIIIBI has to be used. 

One interest of the application of multicomponent in- 
ternal DFT formalism to protons and neutrons is to give 
a fundamental justification to the use of internal density 
functionals for stationary mean-field like calculations of 
nuclei with effective interactions |l2| , showing that there 
exists an ultimate local potential which contains the cm. 
correlations and permits one to reproduce the exact in- 



We have generalized the internal DFT formalism and 
Kohn-Sham scheme to multicomponent self-bound sys- 
tems, treating correctly their fundamental translational 
symmetry. The formalism we obtained applies to the 
description of molecular systems where the nuclei are 
treated explicitly, atomic nuclei and mixtures of 3 He and 
4 He droplets. The main difference with traditional DFT 
is the explicit inclusion of the quantum center-of-mass 



8 Protons and neutrons are both Fermions, so that the two func- 
tionals E x „ should contain exchange and correlations. 

9 3 He are Fermions, so that the corresponding E xc functional 
should contain exchange and correlations. 4 He are Bosons, so 
that the corresponding E x '~ functional should contain only cor- 

relations and the corresponding orbitals <p int should be identical 
(for condensates). 
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correlations (due to the zero-point motion) in the func- 
tional. 

The application to molecular systems where the nuclei 
are treated explicitly permits us to clarify the approxima- 
tions that underly traditional " one kind of particle" DFT 
and KS scheme (i.e. electron-nuclei correlation energy 
neglected) and why this formalism is essentially suited 
for the description of the electrons (only) in the frame 
attached to the cm. of the nuclei. We also set up a 
method to improve the results of traditional DFT in the 
case of molecular systems. 

The application to atomic nuclei and mixtures of 3 He 
and 4 He droplets provides a step towards a fundamental 
justification to the use of effective interactions that are 
functionals of the one-body densities of each species in 
mean-field like calculations. 

Finally, we mention the questions that remain open: 

(i) The search for a general functional that describes 
the cm. correlations and leads to a local KS poten- 
tial is continuing. It would give a numerically ad- 
vantageous alternative to the projection techniques 
used in nuclear physics to restore Galilean invari- 
ance. It could also be used to improve the descrip- 
tion of Helium droplets and nuclei in molecular sys- 
tems. 

(ii) In the case of molecular systems, it appears de- 
sirable to find more satisfactory functional forms 

for the nucleus- nucleus correlation energy E^c [Pint] 
than the simple self-interaction correction to the 
Hartree energy obtained through the crude classi- 
cal (pointlike) approximation. One way could be to 
use semi-classical approximations. It also appears 
desirable to pursue the search of functionals for the 
electron- nuclear correlation energy Eq [Pm > Pint] ■ 

(iii) The generalization this work to the time-dependent 
case is not completely trivial and is under investiga- 
tion, following the considerations of Refs. [32[ and 

m. 

(iv) The question of the non-interacting v- 
representability according to the number of 
orbitals introduced in the KS scheme merits further 
study. 

(v) It would be interesting to apply the same reason- 
ing to rotational invariance to formulate the theory 
in term of the internal density with respect to the 
body-fixed frame (sometimes also called "intrinsic" 
one-body density J46|). Indeed, even if this den- 
sity is not directly observable experimentally, there 
is some indirect experimental evidence for the ex- 
istence of rotational symmetry breaking states in 
nuclear rotational bands. Such a work would cer- 
tainly shed some light on the symmetry breaking 
question. 
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Appendix A: Jacobi coordinates for particles of 
different kinds. 



We consider a system of N particles of different kinds. 
We start from the set of conjugate laboratory coordinates 
{r.;} and {p^}, associated to the masses {rrii}, with i = 
1, ..., N. The Jacobi coordinates are (a = 1, ..., N — 1): 



(i) the cm. position R = r^r Sj=i 



EN 
i=1 mi 



m,rj, conju- 



N 



gated to the total momentum P = X^=i Pi, 

(ii) the relative positions (relative to the cm. of the 
i-1 other particles) £ Q = r a+ i- ^J- . J2t=i TO * r i> 
conjugate to the relative momenta 



r^c + l 

Ei=i m i 



,A 






»=1 P' 



If we note J2i=i mi = M the total mass and 



Pa = 



_ mo, + 1 J27 = i' 



£i=i m i 



the relative mass, we can decompose 



,JV-1 ri 



the kinetic energy as Y,f=i ^~ = JM + ££=i 5jt- 



Appendix B: Simplification of the Jacobi 

coordinates when one kind of particle is much 

heavier than the other. 

We suppose that m' 1 ) >> m^ 2 \ Then, if 7V( 2 ) is not 
much larger than N^ 1 ' and/or R( 2 ' is not very far away 
from R*- 1 ), which is the case for molecular systems, one 
can simplify ([5]) in 



Afd) 



R = — V r (1) = R« 

tt /V(l) 2^ Y 3 "■ ' 



(Bl) 



so that the cm. of the whole system coincides with the 
cm. R^ 1 ' of the particles (1) only. Moreover, Jacobi 
coordinates (HJ) and reduced masses (0 can be simplified 
in 

for ae [1;N (1) -1} : 

k-r^-lf^-tf' andp a = ^-mU, 
a *—* a + 1 

i— i 

for ie [l;iV (2) ] : 
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6v( 



4-l+i 



.(2) 



R« 



'( 2 ) A 

r 4 and \x N m_ x+ i 



,d) 



,(2, 



where the (A^ 1 ' — 1) coordinates £a J appear as Jacobi 
coordinates specific to the particles of kind (1) only, and 

the remaining N^ coordinates r i ' appear naturally as 
the coordinates of each N^ single particles (2) in the 
cm. frame. Thus, we see that we obtain an equiva- 
lent result if we apply Jacobi coordinates to the heavy 
particles only and if we describe each light particle by 
its coordinates in the frame attached to the cm. of the 
heavy particles. As a consequence, the redundant coor- 
dinate problem, and thus the cm. correlations, concerns 
only the heavy particles. This explains why those cor- 
relations should not appear in the exchange-correlation 
energy functional of the light particles (for instance the 
electrons in a molecule). 



Jdr^...d^ldr^...dT%S(R) 



x |V' J :, l t(r 1 " ) ,. 



_W . r ( e ) 



S e ) 



KU)\'fK 



dv 



(e) 



,<<-> 



■ oV ' f(r {e > r 



at(«) ) I 'Pint y 1 1 ' 



(e)/» 



hCW 



M 



r (e) , 



where we have defined the local part of the iV( e )-body 
electronic density as 



dr 



(n) 



dr^MR^) 



(n)^ 



x|^M B) ,- > ra ) ;rW J ...,rW)| 2 .(Cl) 



Appendix C: Exchange-correlation functionals when 
one kind of particle is much heavier than the other. 



This makes it clear that If is no longer explicitly affected 
by (or no longer contains explicitly) the cm. correlations. 
Indeed, the ^(R^™* 1 ) does not directly affect the integral 
over the electronic coordinates. The cm. correlations 



In the molecular case, as R = R( n \ for any function arc included only implicitly through \tjj. 



(e)|2 



(e) 



nates only, we have 



r^ e) ) that depends on the electronic coordi- 



More generally, we consider an operator F that 
acts on the electronic coordinates only (we write 



F(r 



(<0 

i i ■ 



> t nm) m i r } representation). We find 



h = 



(«) 



dr\ 



dv^ 



dv { :l ) dv < t ) ---dv^l ) 8{R)x 



■ N (r, 



</',■ 
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_(") ._(e) 
» 1 AT(")' 1 l ' 



J_(e) „;(e)*/ (e) 

aI AT(«) Vint y 1 ! 1 






■■,'§.>)m 



(e) 
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(e) 
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r (e) V/ e) (r (e) r (e) 1 

1 1 Y< e ) ' Vint \ L l ' ■ • ■ ' 1 nW)i 



U JV< 



r («) ._(e) _(e) < 



(C2) 



where ip^ t is defined such as 



V 



(e)*/_(e) 



dv\ 



(e) 



)F(ri' 



(<■) 



r (e) ^/; (e) fr (e) 



^-•dr^R^*^, 



» .» 



(e) N _ 



■N< 



.)) 



(e) 



,(e) 



)V'™t( r i 



(n) 



» .„(<0 



r («) 



r 



This definition complies with ()Clj) (when F —¥ f). As 
above, we see that J^ is not explicitly affected by the 
cm. correlations (those correlations are included only 

implicitly through the "electronic wave function" f/S)- 
All this makes it clear, from a complementary point 



of view of that presented in appendix [B] that the cm. 
correlations should no longer appear explicitly in the elec- 
tronic exchange-correlation functional (f22j) . Indeed it can 
be rewritten, with the previous notation 
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whose form is similar to that of traditional DFT [l| . Do- ing the same kind of reasoning for the nuclei, we obtain 
I 



E 



(n) 
XC 
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W 1 "' . (n)2 
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i ( ra )*\ 
IViJ ) 



r 



which explicitly contains all the cm. correlations and 
whose form is similar to that of internal DFT ra . 



Appendix D: The classical (point-like) limit in the 
internal DFT formalism. 



To lighten the notation and since it is sufficient for the 
considerations of this paper, we consider only one kind 
of particle in this section. The generalization to many 
kinds of particles is straightforward. 



1. Classical limit for the densities. 

The classical (point-like) limit for a N-body density 
|?/>| 2 is usually obtained by assuming 



W(ri,-..,r,v)| : 



sE" 



N (cl)i 



N\ 



lP W l (r P(i) ), (Dl) 



where the {r^} are the coordinates of the N particles, P 
the possible permutations of those coordinates and 



-,( c 0Vr.\ — 



(r) = 5{v - r. 



(<:V) 



(D2) 



(d) 



are the classical densities of each particle i (the r^ ' be- 
ing their classical positions). Of course, this approxima- 
tion breaks the translational invariance, so that it can be 
fundamentally justified only for densities \ip\ 2 that are 
not translationally invariant (and are symmetric under 
the exchange of two particles, of course). It thus can- 
not be fundamentally justified for the laboratory N-body 
density \ip\ 2 of a self-bound system, nor for its internal 
N-body density |"0mt| 2 : expressed with the particles co- 
ordinates {i"j}, because they arc both translationally in- 
variant. 



For self-bound systems, it is in fact for the cm. frame 
N-body "density" 5(R)\ipi n t(ri, . . . ,tn)\ 2 that the clas- 
sical approximation can fundamentally be justified. In- 
deed, this "density" is obviously not translationally in- 
variant (the <5(R) fixes the cm. in position space and 
amounts to move in the cm. frame) and is symmetric 
under the exchange of two particles (as IV'mtl 2 satisfies 
this symmetry). It is non null only for the {r.;} that sat- 
isfy R = J2i=i r i = (i- c - the {r^} become the cm. 
frame coordinates). 

The same way, for a classical system described in the 
1 p\ c (r P u\) is non null only for 
0. Thus, for self- 



cm. frame, ^rEp n iIiP[ C ° 
the {rj°°} that satisfy £)£=i i[ d) 
bound systems, the proper way to do the classical ap- 
proximation is to replace 

5(K)\^ nt (r u ...,r N )\ 2 -). -Lj2 U liP^ 1 (rp«H D3 ) 



with 



f>j d > = and 



(cl)i 



„(<^ 



Inserting this approximation in the definitions of gTFBJ 
we obtain the classical (point-like) limits of the one and 
two-body internal densities 



N 



pE/M = £>«<')' 



(D5) 



N 



7 $(r. r') = p^(r)^(r') - £ p^ (r)p^f(r'). 
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2. Classical limit for the kinetic energy terms. 

We define the one-body density in {p} space of the 
interacting system as 

Pint(p) = N / dpi ...dp JV (5(P)|V'mt(pi,...,PAr)| 2 (S(p i -p). 



^JV 



where P = J2i=i Pi ^ s nxe d to zero, so that pint is defined 
in the cm. frame. The interacting internal kinetic energy 
can then be written 



JV-l 



v-C r f p 

(lpint\ V TT-l^nt) = / dp — - p int {p). 

'H 2fi a J 2m 

Q — 1 

The one-body density in {p| space of the non-interacting 
(KS) system is defined as [l| 



iV 



pf nt (p) = / , drdrW p(r - r ' )/ft i>l i (r)C t (r'), 

so that the kinetic energy of the non-interacting system 
can be written 

JZi'Pintl-Z-l'Pint) = / d P T~ PintiP)- 



Even if the the one-body densities in the {r} space of 
the interacting and non-interacting systems are the same, 



it is not the case for the one-body densities in the {p} 
space, i.e. Pi n t{p) ¥" PintiP) (the difference is due to the 
quantum correlations). Thus, in the general case, the dif- 
ference between interacting and non-interacting internal 
kinetic energies that appears in Exci see second line of 
Eq. J22J), is not null: 



JV-l 



iV 



(l/>int\ Yl T^HV'mi) - ^2(fint\-^—\<Pint) = 
c=l Zfla i=l Zm 

/2 
dp^{Pint(p)-pL(P))^0- (D6) 



But, since the mean total momenta of the interacting 
and non-interacting systems should be equal, the classical 
(point-like) limit (obtained supposing that the system is 
perfectly localized in both {r} and {p} spaces) gives 



Pint(p) = PintiP) 



M) 



N 

£ 

t=i 



6(p 



Mh 



where the p^ ' are the classical momenta of the particles 

in the cm. frame, which satisfy Ylt=i Pi = ( m ^ ne 
stationary case, pf = of course). Thus, in the classi- 
cal limit, the "kinetic energy exchage-corrclations" part 
of Exc (i- G - the second line of Eq. (|2"2"j) ). which contains 
explicitly the cm. correlations in {r} representation, dis- 
appears. 
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